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Abstract 

We report the first results from tlie Clusters Around Radio-Loud AGN (CARLA) program, a Cycle 
7 and 8 Spitzer Space Telescope snapshot program to investigate the environments of a large sample 
of obscured and unobscured luminous radio-loud AGN at 1.2 < z < 3.2. These data, obtained for 
387 fields, reach 3.6 and 4.5 /im depths of [3.6]ab = 22.6 and [4.5]ab = 22.9 at the 95% completeness 
level, which is two to three times fainter than L* in this redshift range. By using the color cut [3.6]- 
[4.5] > —0.1 (AB), which efficiently selects high-redshift {z > 1.3) galaxies of all types, we identify 
galaxy cluster member candidates in the fields of the radio-loud AGN. The local density of these 
IRAC-selected sources is compared to the density of similarly selected sources in blank fields. We find 
that 92% of the radio-loud AGN reside in environments richer than average. The majority (55%) of 
the radio-loud AGN fields are found to be overdense at a > 2 cr level; 10% are overdense at a ^ 5 cr 
level. A clear rise in surface density of IRAC-selected sources towards the position of the radio-loud 
AGN strongly supports an association of the majority of the IRAC-selected sources with the radio- 
loud AGN. Our results provide solid statistical evidence that radio-loud AGN are likely beacons for 
finding high-redshift galaxy (proto-) clusters. We investigate how environment depends on AGN type 
(unobscured radio- loud quasars vs. obscured radio galaxies), radio luminosity and redshift, finding 
no correlation with cither AGN type or radio luminosity. We find a decrease in density with redshift. 
consistent with galaxy evolution for this uniform, flux-limited survey. These results are consistent with 
expectations from the orientation-driven AGN unification model, at least for the high radio luminosity 
regimes considered in this sample. 

Subject headings: galaxies: active — galaxies: clusters: general — galaxies: high-redshift — infrared: 
galaxies — techniques: photometric 



L INTRODUCTION 

The Infrared Array Camera ('IRAC: lFai"io et al.ll200l) 
of the Spitzer Space Telescope is an incredibly sensitive 
tool for finding and studying massive galaxies at high 
redshift. For stellar populations formed at high redshift, 
negative fc-corrections provide a nearly constant 4.5/im 
flux density over a wide redshift range — e.g., an L* 
galaxy formed at Zf = 3 will have [4.5] ^ 17 at 0.7 ^ z ^ 
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2.5, which is sufficiently bright that it is robustly seen in 
even 90 sec integrations with Spitzer. Several teams have 
been exploiting this capability to identif y large samples 
of galaxy clusters at z ^ 1 (e.g., [E iscnh ardt et al.l [20081: 
Wflsoii et al.ll2009l:lPapovich et al.ii2010i:iGalametz et al.l 



2010aHMuzzin et al.ll2Q13| ). TheSpitzer mid-infrared se- 
lection is capable of flnding both the very rare, massiv e 
evolved clusters at z ^ 1.75 (e.g.. [Stanford et ani2012| ). 
as well as the more numerous gro ups or forming clus- 
ters at similarly h igh redshifts (e.g.. lZeimann et al.ll2012l : 
iGobat et al.ll201l[) . Teams are now using mid-IR-selected 
clusters for a range of important studies such as using 
the colors, sizes, and mid-IR cluster galaxy luminosity 
function to pr obe the formation epoch of niassive cluster 
galaxies fe.g. iMancone et al.l 120101 120121: JRettura et al " 



2011 



: iLidman et al.ll2012t iMei et alj|20ia ISnvder eFa! 



2012f ). measuring the galaxy cluster co rrelation function 
out to z '^ 1.5 (jBrodwin et al.l 120071) . probing e volu- 
tion in the a — T^ correlation ( Brodwin et al.l[2011[ ). and 
targeted cosmological studies of z > 1 Type la. SNe 
in dust-free cluste r environments ([Dawson et al.l [20091 : 
ISuzuki et al"][20Tl) . 

Much of the work to date on high-redshift galaxy 
clusters comes from field surveys, which provide both 
a strength and a weakness. Uniformly selected galaxy 
cluster samples have the power to probe basic cosmo- 
logical parameters by measuring the growth of struc- 
ture. However field s urveys — both in the mid-IR (e.g., 
[Gettings et al.l [2012[ ) and at other wavelengths (e.g., 
Sunyaev-Zeldovich and X-ray survevs: .Vanderlinde et al.l 
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[Mot iFassbender eFall [20lTt iReichardt efaLl [20T2I ) — 
find few clusters at tlie fiigliest redsliifts, z ^ 1.5. Wliile 
measuring tlie cluster space density is required to use 
these systems as cosmological probes, many key galaxy 
cluster studies d o not require such kn o wledge (e.g., 
Krick et all 120081: IGalametz et al.l I2009t ISuzuki et all 



2012t iMarthii et al.(l2013[ r 



In order to efficiently identify the richest environ- 
ments at yet higher redshifts, targeted searches for high- 
redshift galaxy clusters have many advantages. To- 
wards that goal, literature that stretches back nearly 
50 years shows that powerful radio-selected AGN pref- 
erentially reside in luminous red sequence galaxies (e.g. , 
Matthews! [I96l iBest et al.lll998t IVenemans et all 120021: 



Hickox et all 120091 [Griffith fc SternI 120*^ ? Indeed, tar- 
geted searches for high-redshift clusters and proto- 
clusters around powerful high-redshift radio galaxies 
(HzRGs) have proven very successfu l and have a rich 
literature (e.g.. IPentericci et al.ll2000l: iStern et al.ll2003l: 
Kurk et all 120041: IVenemans et al.l 120071 : iDohertv et al.l 



Kurit et ai.l i:^UU4i: ivenemans et ai.l \i\i\iix. luonerty et ai.l 
20101: IGalametz et al.l 120091 I2010bl: iHatch et al.l l20Tl[ 



Matsuda et al.ll201lHMavo et al.l | 2012D 



In a pilot study by our team, IGalametz et al.l (J2012D 
used a counts-in-ccll analysis to identify overdensities 
of IRAC-sclectcd high-redshift galaxy candidates in the 
fields of 72 HzRGs from the Spitzer High-Redshift 
Radio Galaxy p r ogr am (SHzRG; iSevmour et al.l 120071 : 
IDe Breuck etalj [201 Oi ). The HzRGs in that study 
spanned a range of 1 .2 < z < 3. Using relat ively shallow, 
- 120 s IRAC data. IGalametz et all (|2012D showed that 
radio galaxies preferentially reside in medium to dense 
regions, with 73% of the targeted fields denser than av- 
erage. Apart from six newly discovered cluster candi- 
dates, several known (proto-)clusters were recovered as 
overdense fields in that analysis, and we have recently 
spectroscopically confirmed one of the more promising 
candidate fields as a new high-redshift cluster at z = 2.02 
(Galametz et al., in prep.). 

In this paper, we describe the first results from our 
^ 400 hr Warm Spitzer snapshot program called Clus- 
ters Around Radio-Loud AGN, or CARLA[3, a Cycle 7 
and Cycle 8 program which takes advantage of Spitzer''s 
impressive sensitivity to identify massive galaxies at high 
redshift. To date, CARLA has targeted nearly 400 radio- 
loud AGN, including 187 radio-loud quasars (RLQs) and 
200 HzRGs. Targets were uniformly selected in radio lu- 
minosity over the redshift range 1.3 < z < 3.2. Simila r 
to the successful pilot study bv IGalametz et al] (J2012D . 
CARLA allows, for the first time, a sensitive, systematic 
study of the environments of a large sample of powerful 
radio-loud sources over a wide redshift range. We iso- 
late high-redshift galaxy candidates using IRAC colors 
and quantify the environments of the targeted AGN as 
compared to similar-depth blank-field surveys. CARLA 
also allows us to investigate how environment depends on 
radio-loud AGN properties such as redshift, radio lumi- 
nosity and AGN class (e.g., type-1, or unobscured RLQs 
vs. type-2, or obscured HzRGs). A companion paper, 
Wylezalek et al. (in prep.), investigates the luminos- 
ity function of CARLA-selectcd candidate galaxy cluster 

i"* As of March 2013, CARLA targets arc still being scheduled; 
this manuscript details all observations obtained through February 
2013. 



40 



"20 




29.5 



29.0 



X 

g 

-i: 28.5 




28.0 



27.5 
1 









+ 

o 














20 40 









I— 1 


II 



1.5 



2.0 



2.5 



3.0 



3.5 



50 



-50 









+ 






+ 


ii*r|*- 


v/ V: 


* f+ + 


* * V 




+ ++* + 



HzRG 
RLQo 



+ 
"+ ++ 1 



100 200 

RA (J2000) 



300 



Fig. 1. — Top panel: Rest-frame 500 MHz luminosity vs. red- 
shift. As shown by their distributions, both subsamples (HzRGs 
and RLQs) sample the targeted LsqomHz ~ ^ plane in a relatively 
uniform manner. Bottom panel: Distribution of the CARLA radio- 
loud AGN across the sky. Most of the RLQs are chosen from the 
SDSS and are thus mainly located in the Northern hemisphere, 
while the HzRGs cover the full extragalactic sky. In both panels, 
green crosses show the HzRGs in our sample, while orange dia- 
monds show the RLQs. 



members. 

The paper is organized as follows: § 2 describes how 
the sample of HzRGs and RLQs was selected and § 3 
gives details on the observations, data reduction, catalog 
generation and the completeness limit of the survey. We 
describe the analysis of the environments of the CARLA 
targets with respect to galaxy type, redshift and radio 
luminosity in § 5 and draw conclusions from the results 
in § 6. Section 7 summarizes our results. Throughout 
the paper wc assume Hq = 70 km s~^ Mpc~^, fim — 0.3, 
Q\ = 0.7. All magnitudes and colors arc expressed in 
the AB photometric system unless stated otherwise. 

2. THE CARLA SAMPLE 



High-redshift galaxy clusters around radio-loud AGN 



A key goal of the CARLA survey is to investigate 
how environment depends on AGN classification, and 
we therefore have targeted a matched sample of unob- 
scured (e.g., type-1) and obscured (e.g., type-2) radio- 
loud AGN. Type-1 AGN, specifically broad-lined quasars 
or RLQs for the high-luminosity, radio-loud sources tar- 
geted here, show signatures of radiation originating from 
very close to the central black hole. In contrast, the cen- 
tral black hole is obscured by dense absorbing material 
in type-2 AGN, or HzRGs for the high-luminosity, high- 
redshift, radio-loud sample obs erved here. In the stan- 
dard AGN unification model (jUrrv fc PadovamI I1995D 
this obscuration is produced by an optically thick cir- 
cumnuclear torus and differences in optical appearance 
only depend on the orientation angle at which the AGN is 
observed, not due to intrinsic differences between the two 
AGN types. A clear prediction is that both AGN types 
should reside in identical environments. Recently, an al- 
ternative paradigm has been gaining attention in which 
type-2 AGN evolve into type-1 AGN as their merger- 
triggered, dust-obscured AGN be come more powerful 
and clear out the environment (e.g. iHopkins et al.|[2006l 
and references therein). This dynamical model would 
mean that type-1 and type-2 AGN trace different phases 
in the evolution of powerful AGN, potentially associated 
with different environments. On the other hand, if the 
timescale associated with the AGN breaking out of its 
obscuring cocoon is much shorter than the timescale on 
which the environment develops, then no dependence of 
environment on AGN type would be expected in the dy- 
namical model. 

The CARLA sample consists of 387 radio-loud AGN, 
including 187 RLQs and 200 HzRGs (sample size as of 
2013 February 20). The matched RLQ and HzRG sam- 
ples have been deliberately designed to be of similar size 
with similar redshift and radio luminosity ranges. HzRGs 
are defined as powerful, high-redshift radio sources with 
rest-frame 500 MHz radio luminosities isooMHz ^ 10^^'^ 
W Hz~^ and have optical spectra free of broad emis- 
sion lines. This classification as narrow-line, i.e. type- 
2 AGN is based on the optical spectroscopy used to 
determine their redshifts. The H zRGs were selected 
from the updated compendium of iMilev fc De Breuc3 
()2008() extended down to z = 1.3 using both flux- 
limited radio surveys (e.g. M RC, 3C, 6C, 7C) and 
ultra-steep spect r um su rveys (e.g. lRottgering et al.lll997l : 
iDe Breuck et al.l l2001| ) . The sample contains all 146 
powerful HzRGs known as of 2010 at 2 < z < 3.2. 

To calc ulate Jvbqqmhzi we used the procedure de- 
scribed bv lMilev fc De Breuckl ()2008i ) . In short, we cross- 
correlated the Hz RGs with the 1 4 GH z NRAO VLA Sky 
Survey (NVSS; ICondon et all [1991 and the 74 MHz 
VLA Low-Frequency Sky Survey (VLSS; iCohen et al.l 
|2007| ). using a 30" search radius. These two radio sur- 
veys cover a very large part of the extragalactic sky, and 
have well-matched spatial resolution, thereby avoiding 
concerns with missing fiux. The rest-frame 500 MHz fre- 
quency was deliberately chosen to ensure that the fluxes 
were interpolated rather than extrapolated between the 
two radio surveys. If no VLSS detection was avail- 
able, we assumed a. spect ral index a = —1.1 (Si, oc i/", 
IMilev fc De Breuckl [20081) . 

To create the RLQ sample, we selected all optically 
bright {Mb < —26.5) quasars in the Sloan Digitized Sky 



Survey (SDSS; ISchneider et^[20lol) and the 2 dF QSO 
Redshift Survey f2dFQZ: iCroom et al.l [200l . where 
the spectroscopy-based classifications were taken from 
the original catalogues. For the SDSS quasars the AIb 
were taken from the original catalogue, and for the 2df 
quasars we calculated Mb from the original catalogs 
assuming a quasar SED represented by a power law 
(/y oc :y"«^Q , aB.LQ = 0.5 ~ as used for the SDSS quasars) 
and we corrected for Galactic dust extinction. 

To identify RLQs, we then correlated this sample with 
the NVSS using a 10" search radius. As only '-20% of the 
RLQs are detected in the VLSS, we could not calculate 
the radio spectral indices in a systematic way. Moreover, 
RLQs may contain time- variable, Doppler beamed emis- 
sion, which would render the spectral indices from sur- 
veys observed several years apart unreliable. Wc there- 
fore assumed a fixed a = —0.7 for all RLQs to calculate 
their Lsoomhz- We tested and justified this choice by 
computing the mean and median value for the 133 out of 
717 RLQs where the VLSS-NVSS spectral indices could 
be derived, and find —0.65 and —0.74, respectively. On 
average, the radio luminosities calculated with a fixed a 
only differ from the precise radio luminosities by 3% and 
this scatter will not have any influence on our results and 
conclusions. 

The final RLQ sample was obtained by considering 
only sources in regions with low Galactic extinction 
{E{B — V) ^ 0.1) and by matching to the redshift and 
radio power distribution of the HzRG sample. The up- 
per panel in Fig. [T]shows rest-frame Lsoomhz vs. redshift 
while the lower panel shows the distribution of our tar- 
gets across the sky. This sample selection allows us to 
study evolutionary trends as well as trends with radio 
luminosity for both types of AGN without any biases 
beyond the selection being restricted to the most radio 
powerful sources. 

3. DATA 

3.1. Observations 

CARLA observed the lower redshift (1.3 < z < 2) 
sources with the Spitzer IRAC camera at 3.6 and 4.5 
/xm (referred to as IRACl and IRAC2) during Cycle 7 
with total exposure times of 800s in IRACI and 2000s in 
IRAC2. Slightly deeper observations were obtained for 
the higher redshift (2 < z < 3.2) portion of the sam- 
ple during Cycle 8 with total exposure times of 1000s 
in IRACl and 2100s in IRAC2. IRAC employs 256 x 
256 InSb detector arrays with 1'.'22 pixels to map out 
a 5.2' X 5.2' field. The IRACl and IRAC2 data were 
obtained simultaneously over slightly offset fields, allow- 
ing CARLA to tailor the exposure times in each band 
for our intended science. Observations were all obtained 
with dithered 100 s observations. In order to save the 
^ 5 min overhead per target had the differing integra- 
tion times for IRACl and IRAC2 been split across two 
Astronomical Observation Requests (AORs), we instead 
implemented this observing strategy in a single AOR us- 
ing the the "Fixed Cluster Offset" mode in array coor- 
dinates. In AB magnitudes, the total exposure times 
give similar depths for both channels (see section 3.4) 
with the IRAC2 observations being slightly deeper. The 
total exposure times, which arc slightly longer for the 
high-redshift portion of the sample (observed during Cy- 
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clc 8), reach depths two magnitudes below m* over the 
whole redshift range (see Section 3.4). This allows us to 
robustly measure the faint-end slope of the cluster lumi- 
nosity function out to high redshift, which is the focus of 
our second CARLA paper (Wylezalek et al., in prep.). 

3.2. Data Reduction 

The basic calibrated data were reduced and mosaicked 
using the MOPEX package (Makovoz & Khan 200£i) and 
resampled to a pixel scale of 0'.'61. The MOPEX outlier 
(e.g., cosmic ray, bad pixel) rejection was optimized for 
the regions of deepest coverage in the center of the maps, 
corresponding to the location of the targeted AGNs. Due 
to artifacts and bright foreground stars, the maximum 
coverage is not reached in every part of the field. Regions 
with ^ 85% of the maximum coverage were ignored for 
source extraction and further analysis. 

3.3. Source Extraction 

Source extraction wa s performed using SExtractor 
(jBertin fc Arnoutill996D in dual image mode using the 
4.5 /xm frame as the detection image. W e used the 
IRAC -optimized SExtractor parameters from lLacv et al.l 
()2005() . We measured flux densities in 4" diameter aper- 
tures, converting from the native MJy stcr~^ units of 
the images to fiJy pixel" ^ by multiplying with the con- 
version factor of 8.4615 /xJy pixel"^ / (MJy ster~^) for 
our 0'.'61 pixel scale. We derived empirical aperture cor- 
rections from the curve of growth of bright, isolated stars 
in the CARLA survey. These multiplication factors, 1.42 
for IRACl and 1.45 for IRAC2, correct the photometry 
to 24"-radius apertures. These values are systematically 
~ 17% larger than the aperture corrections listed in the 
Spitzer Instrument Handbook but are consistent w ith the 
aperture corrections found bv lAshbv et al.l (|2OO90 for the 
Spitzer Deep, Wide-Field Survey. The discrepancy with 
the Instrument Handbook is most likely due to differ- 
ences in the backgroun d determination an d slight differ- 
ences in aperture radh (|Ashbv et al.ll2009l ). 

3.4. Completeness Limit 

We determined the 95% completeness level by compar- 
ing the CARLA number counts to number counts from 
the Spitzer UKIDSS Ultra Deep Survey (SpUDS, PI: J. 
Dunlop), a deep Cycle 4 Spitzer Legacy program cover- 
ing 1 deg^ in the UKIDSS UDS field with IRAC and the 
Multiband Imagin g Spectrometer aboard Spitzer (MIPS; 
iRieke et al.ll2004 ). SpUDS reaches greater sensitivities 
than CARLA, with IRACl and IRAC2 3 a depths of 
^ 1 /xJy (magAB — 24). Although SpUDS catalogs are 
publicly available, we performed source extraction on the 
final mosaics in the same way as for CARLA to ensure 
uniformity in our analyses. In particular, the public cata- 
logs are cut at a depth of 8 a, i.e. shallower than CARLA, 
and thus would not allow us to accurately measure the 
completeness of the CARLA catalogs. The aperture cor- 
rections we derive for bright, isolated sources in SpUDS 
are consistent with the ones derived for the CARLA sur- 
vey. We created IRACl and IRAC2 SExtractor catalogs 
using single image mode. The differential number counts 
are shown in Fig. [H As planned, the IRACl and IRAC2 
observations match very well in depth. We determine 
a 95% completeness limit of 3.45 and 2.55 fj,Jy relative 



to the SpUDS survey, corresponding to limiting magni- 
tudes of [3.6] = 22.6 and [4.5] = 22.9 for the observations. 
These values are adopted as flux density cuts for all sub- 
sequent analyses. 

4. THE ENVIRONMENTS OF RADIO-LOUD AGN 

4.1. IRAG-Selected Sources 

Wc use a counts-in-cell analysis of color-selected 
IRAC sou r ces to identify high-redshift overdcnse fields. 
iPapovichI (|2008[ ) showed that an IRAC color cut of 
[3.6] — [4.5] > —0.1 efficiently isolates galaxies regard- 
less of age and galaxy type at z > 1.3. This crite- 
rion makes use of the 1.6 fxin bump, a prominent spec- 
tral feature in the spectral energy distributions of galax- 
ies that does not depend on the evolutionary stage of 
the galaxy. The bump is caused by a minimum in the 
opacity of the H^ io n which is present in the atmo- 
spheres of cool stars (jJohnI 1198 8') and is a commonly 
used i ndicator of redshift (e. g., Simpson & Eiscnhard^ 
119991: iSorbafc Sawickill2010D . The 1.6 fim bump en- 
ters the IRAC bands at z ~ 1 and provides a redder 
[3.6] — [4.5] color for high-redshift galaxies as compared 
to low-redshift (z < 1) galaxies (see Figl3]). Contam- 
ination can come from strongly s tar-forming galaxies 
at 0.2 < z < 0.5 (jPapovichll2008D . cool brown dwarfs 
i ■ 

( 
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200 



200. 



and powerful AGN at all redshifts 
However, these contaminations are 
not expected to be signific ant (for a detailed discussion 
see iGalametz et al.ll2012D . More than 90% of objects 
satisfying this simple color criterion will be at z > 1.3 
(jPapovich 200^ . In the following we refer to these high- 
redshift, mid-IR color-selected sources simply as TRAC- 
selectcd sources'. 

4.2. Comparison to Blank Fields 

We use the SpUDS survey to derive a mean blank field 
density of IRAC-selected sources. The catalogs have 
been cut at the CARLA depth, i.e. all sources fainter 
than the CARLA 95% completeness limit are neglected. 
We then placed 436 non-overlapping circular cells of 1' 
radius onto the field and measured the density of IRAC- 
selected sources in these cells. None of the CARLA tar- 
gets lie in the SpUDS survey. 

We then measured the number of IRAC-selected 
sources in an aperture of 1' radius centered on the radio- 
loud AGN. At 1 < z < 3, 1' corresponds to -- 500 
kpc, matching typical cluster sizes for z > 1.3 mid- 
IR selected clusters with log(M2oo/Af0) ^ 14.2 (e.g. 
lBrodwinet"aII 120111) . The surface density of IRAC- 
selected sources in the CARLA fields takes into account 
the low coverage areas (see Section 3.2). Table [1] lists 
the names, positions, redshifts, radio luminosities and 
surface density of IRAC-selected sources in the CARLA 
fields sorted by right ascension. The full table is avail- 
able online. The upper panel in Fig. 2] shows the dis- 
tribution of surface densities for IRAC-selected sources 
in the SpuDS and CARLA fields. A perfect, structure- 
free blank field would show a Gaussian distribution of 
densities peaking at the mean blank field density. As ex- 
pected, both the SpUDS and CARLA distributions are 
asymmetric with a high-density tail since even the '^ 1 
deg^ SpuDS field contains large scale structure. For ex- 
ample, a high-density tail of IRAC-selected sources is also 
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Fig. 2. — Number counts for the SpUDS and CARLA surveys for sources detected at 3.6 fira (left) and 4.5 fim (right). These analyses 
rely on catalogs independently built for detecting sources in each band. The 95% completeness limit in each band is derived by comparing 
the CARLA number counts to those from the SpUDS survey, which gives completeness limits of 3.45 fi3y for IRACl and 2.55 /ijy for 
IRAC2, as shown by the vertical dotted lines. 
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Fig. 3. — Illustra t ion of the IRAC color criterion. Left: Predicted evolution of IRAC2 magnitude as a function of redshift. We plot the 
IBruzual fc CharloU 1 12003 ') models for different formation redshifts Zf assuming a Salpeter initial mass function and a single exponentially 
decaying burst of sta r formation with r = 0.1 Gyr, normalized to the observed m* for galaxy clusters at 2 = 0.7, [4.5] = 16.75 (Vega; 
IMancone et al.ll201CII ). The models were generated with the model calculator EzGal l IMancone fc Gonzalej[2012l ). A negative fc-correction 
leads to an almost flat evolution of the IRAC2 magnitude at 2 ^ 0.7. Right: The shift of the 1.6 ^m bump across the IRAC bands leads 
to a red IRAC color for galaxies at 2 > 1.3. 



TABLE 1 
The CARLA sample (in R.A. order) 



Name 



R.A. 
(J2000) 



Dec. 
(J2000) 



log(L500MHz) 2 

(W Hz"-'-) (arcmin^^) 



AGN Type 



U9S0003-019 


00:06:11.00 


-01:41:50.2 


1.54 


27.86 


10.5 


HzRG 


PKS.0011-023 


00:14:25.00 


-02:05:56.0 


2.08 


28.15 


10.2 


HzRG 


MRC.0015-229 


00:17:58.20 


-22:38:03.8 


2.01 


28.32 


15.6 


HzRG 


BRL0016-129 


00:18:51.40 


-12:42:34.6 


1.59 


29.02 


15.3 


HzRG 


MG0018+0940 


00:18:55.20 


-1-09:40:06.9 


1.59 


28.39 


13.4 


HzRG 



This table is published in its entirety in the electronic edition of ApJ; a portion is shown here for guidance regarding its form and content. 
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Fig. 4. — Upper panel: Histogram of the densities of IRAC- 
selected sources in the CARLA fields and the SpUDS survey. The 
Gaussian fit to the low-density half of the SpUDS density distribu- 
tion is shown by the dashed black curve, giving SgpUDS = 8.3±1.6 
arcmin"^. The vertical dotted line corresponds to EspUDS + Sc 
density. We find that 37.0% of the CARLA fields are over- 
dense compared to this surface density while only 8.6% of the 
SpUDS fields have surface densities of IRAC-selected sources above 
that value. Lower panel: Cumulative density distribution of the 
CARLA and SpUDS fields. Only 18.7% of the SpUDS fields 
arc denser than SspUDS + 2(TSpUDS in contrast to 55.3% for the 
CARLA fields. Only 0.7% of the SpUDS fields are denser than 
SspUDS + So'SpUDS while 9.6% of the CARLA fields are denser 
than this cut. 



seen in the Lockman H ole field iP apovich 200 ^. Follow- 
ing th e met hodology oflPapovichl (,2008. ) , .Galametz et al.l 
(|2012D and IMavo et al.l (J2012D . we fit a Gaussian to the 
low-density half of the distribution (i.e. the distribu- 
tion of the lower density regions determined with re- 
spect to the maximum of the distribution) to determine 
the peak density SspUDS = 8.3 arcmin~^ and the stan- 
dard deviation of the low-density half, crspUDS =1-6 
arcmin-2. We find that 92.0% of the CARLA fields arc 
denser than the SpUDS peak density, 55.3% arc ^ 2 ct 
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Fig. 5.— Postage stamps of the fields of MRC1217-276 and 
J105231.82-f08060, a radio galaxy and a radio-loud quasar at 
z = 1.90 and z = 1.64, respectively. The IRAC color-selected 
sources within a 1 arcmin radius, shown with the large dashed cir- 
cle, are marked with small solid red circles. The radio-loud AGN 
are marked by the orange squares. 

overdense and 37.0% are ^ 3 ct overdense. At the 5 
a level, 9.6% of the CARLA fields arc still overdense. 
The corresponding values for the SpUDS field arc signif- 
icantly lower, with 18.7%, 8.6%, and 0.7% of the fields 
being denser than SspUDS + SctspUDS, SspUDS+3aspUDS, 
SspUDS + ScspUDS, respectively (see Fig. |4l). The densest 
CARLA field, a HzRG at z = 2.00, shows a density of 
almost SspUDS + lOcr. A two-sided Kolmogorov-Smirnov 
(K-S) test gives a < lO'^*^ probability that the CARLA 
and SpUDS surface density distributions are drawn from 
the same underlying distribution. In addition, a Mann- 
Whitney U test rules out that the two distributions have 
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the same mean. 

As the CARLA fields are distributed over the whole 
sky, we compute the Spearman rank correlation coeffi- 
cient p to test the dependence of density on Galactic 
latitude to rule out any bias due to stellar contamina- 
tion. We find p = 0.01, implying a 21% probability of 
a correlation. According to standard interpretations of 
the Spearman rank correlation test, this is considered as 
e vidence for the lack of a correlation. 

iGalametz et al.l (|2012( ) also used a counts-in-cell anal- 
ysis to identify overdensities of IRAC-selected sources 
among 72 HzRGs in a redshift range of 1.2 < z < 3, 
albeit with observations two magnitudes shallower than 
CARLA. They found that 73 ± 12% of the fields were 
denser than the mean blank-field density, derived from 
the SWIRE survey They also found that 23 ±7 % of 
the fields were ^ 2 cr ove rdense. We cut the C ARLA 
survey at the depth of the IGalametz et al.l ()2012l ) study 
to check for consistency with this previous work. The 
SpUDS mean blank-field density, derived from a Gaus- 
sian fit to the low density half of the distribution, gives 
2.1 ± 0.7 IRAC-selected sources per arcmin^, in good 
agreement with the bl ank field mean density derived by 
IGalametz et al.l (|2012[ ). Sswire = 2.8 ± 1.0 arcmin^^ 
Most of the CARLA fields are denser than the blank 
field, with 79 ± 9% denser than the SWIRE surface den- 
sity, consistent with IGalametz et al.l (|2012[ ). At the high 
end of the density distribution we find more fields, with 
51 ± 7% being ^ 2 a overdense. This may, however, be 
due to small number statistics in IGalametz et aP ()2012l ) 

In Fig. 2] we show the positi on of a known cluster , 
GIG 0218.3-0510 at z = 1.62 (IPapovich et al.l l2010l) . 
in the density distribution. This cluster lies in the 
SpUDS survey and is centered on one of our random 
apertures placed onto the survey. The aperture cover- 
ing thi s cluster is the denses t one in the SpUDS sur- 
vey. In IGalametz et al.l (|2012f ) , the IRAC color criterion 
has been successfully tested by recovering known high- 
redshift (proto-) clusters such as PKS1138— 262 {z = 
2.16, the Spiderwcb galaxy). In addition, we recently 
confirmed a cluster in the field of MRC0156-252, a radio 
galaxy at z = 2.02, which is 2 ct overdense in the CARLA 
data (Galametz et al., in prep.). 

These results confirm that the color criterion is very 
effective at identifying high-redshift clusters and that 
the densities found in this work are sufficient to qual- 
ify many of the fields as promising (proto-)cluster can- 
didates. To demonstrate the quality of our data. Fig. 
[5] shows the IRAC2 images of two CARLA sources, the 
quasar J105231.82-f08060 {z = 1.64; 6.8 a overdense) 
and the radio galaxy MRC1217-276 {z == 1.90; 5.3 a over- 
dense). IRAC-selected sources within 1' of the radio- loud 
AGN are indicated, illustrating the diversity of environ- 
ments. While the former shows a central concentration 
of IRAC-selected sources close to the radio-loud AGN, 
the latter shows a more diffuse overdensity, perhaps sug- 
gestive of an earlier stage in the cluster formation. 

4.3. Radial Distribution of IRAC-Selected Sources 

We investigate how the surface density of the IRAC- 
selected sources depends on the distance to the radio- 
loud AGN. In Fig. [H] we show a stack of the surface 
density of IRAC-selected sources in annuli around the 
radio-loud AGN for both the CARLA and the SpUDS 
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Fig. 6. — Surface density of the IRAC-selected sources as a func- 
tion of distance to the radio-loud AGN for the CARLA survey. 
We show a stack of all CARLA fields that are denser and less 
dense than SspUDS (black asterixes and red crosses, respectively) 
and a stack of fields denser than SgpUDS "I" ScrgpUDS (blue dia- 
monds). The radial distribution depends on the overall density 
of the CARLA fields with a sharp rise towards the center for the 
densest CARLA fields and an almost homogenous distribution for 
the lowest density fields. The grey dashed line corresponds to the 
SpUDS peak density and is slightly lower than the SpUDS radial 
distribution which is computed from a stack of all SpUDS fields, 
including both the low- and the high-density fields. 

survey. As expected, the radial distribution of IRAC- 
selected sources in the SpUDS survey is homogenous and 
is equivalent to a median blank field surface density. It 
should be noted that this median is dependent on large- 
scale structure and existing clusters in SpUDS and is 
therefore slightly higher than the typical blank field den- 
sity derived from the Gaussian fitting to the lower part 
of the density distribution (see Fig. 2]). For the dense 
CARLA fields, the surface density peaks towards the po- 
sition of the radio-loud AGN in the center of the image. 
This rise is even more prominent for the CARLA fields 
overdense at the 3 <t (E = SspUDS + SfXspUDs) level than 
for the CARLA fields that are simply denser than the 
peak SpUDS density. This means that, on average, the 
denser a field the more IRAC-selected sources are found 
very close to the radio-loud AGN. This result clearly 
shows that for the dense CARLA fields, the radio-loud 
AGN is a good beacon for the center of the overdensity 
of IRAC-selected sources. 

We also test that this behavior is not simply due to 
the CARLA sample centering apertures on known mas- 
sive, high-redshift galaxies while the SpUDS compari- 
son sample targeted random positions in that field. We 
identify a subsample of SpUDS sources with similar mid- 
infared colors and magnitude distributions to the HzRGs 
in CARLA; the CARLA RLQs are too dominated by 
their AGN in the mid-infrared to use for this compari- 
son analysis. We find that the density of IRAC-selected 
sources for this matched sample does not show a promi- 
nent peak towards the center, and is consistent with the 
SpUDS blank field density at all radii, thus further sup- 
porting our primary result that powerful radio-loud AGN 
are efficient beacons for high-redshift structure. 

Finally, we also investigate the radial density distribu- 
tion for all CARLA fields that are less dense than the 
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SpUDS peak density to rule out systematics in our anal- 
ysis. The rise towards the center fully disappears and is 
consistent with the SpUDS peak density shown by the 
grey dashed line. It is important to note that even at 
largest radii probed, the CARLA fields show an excess 
of IRAC-selected sources suggesting that the extent of 
the o verdensity is larg e r than the area covered in our im- 
ages. iVenemans et al.l ()2007() measures the size of (proto- 
) clusters by analyzing the radial extent of Lya emitters 
in the fields of radio galaxies at z ~ 3. This work showed 
that only at a ~ 1.75 Mpc (~ 3.8 arcmin) distance from 
the radio galaxy does the density of Lya emitters in the 
field become consistent with the field density. They con- 
clude that the sizes of (proto-) clusters as inferred from 
the distribution of Lya emitters at z ~ 3 is roughly 2 
Mpc. Due to the limited size of the CARLA fields we 
cannot measure the overall extent of the overdensity but 
we show that it extends beyond 2 arcmin {'^ 1 Mpc for 
the redshifts probed). The clear rise towards the center 
at distances < 60 arcsec, however, justifies our choice of 
1 arcmin apertures for the counts-in-cell analysis. 

4.4. Dependence on AGN Type 

We compare the fields around the HzRGs and the 
RLQs separately to investigate any difference in their 
environments. Fig. [7] shows the density distributions 
for the HzRGs and the RLQs in our sample. We pro- 
ceed in the same way as earlier and fit a Gaussian to 
the low density part of the histograms to find the mean 
and the standard deviation of the distributions. The 
best fit is found for a mean surface density of IRAC- 
selected sources of 11. 6 ±2. 2 arcmin"^ for the HzRGs and 
10.8±1.7 arcmin"^ for the RLQs. The fitted means agree 
within the uncertainties. A two-sided K-S test shows the 
two distributions to be statistically similar, with only a 
38% chance of the two distributions being different. The 
Mann- Whitney U test only gives a 23% chance of the two 
means of the distributions being different. We repeat the 
analysis for various redshift bins to test if this result is 
only true for specific redshifts, finding similar results for 
all tested redshift bins. We discuss the implications of 
environment not correlating with AGN obscuration in 
Section 5. 

4.5. Dependence on Redshift 

Figure |8] shows the environments of radio-loud AGN 
as a function of redshift. Spearman rank correlation co- 
efficients suggest a dependence on redshift with a prob- 
ability of 99.99%. This correlation is mainly driven by 
fewer data points at z < 2.0 and E < 10 arcmin"^ and 
at z > 2.5 and S > 12 arcmin"^. We therefore also test 
for a correlation in these two redshift bins (1.3 < z < 2.5 
and 2.5 < z < 3.2); the Spearman rank correlation co- 
efficients suggest a correlation with a probability of 99% 
and 82%, respectively and confirm that the dependence 
of density with redshift is true for all redshifts probed 
here. The density of the environments around the RLQs 
seems to be more dependent on redshift than those of 
the HzRGs with probabilities of 99.0% and 83.0% of a 
correlation, respectively. However, this could be due to 
slightly different sampling at intermediate redshifts. If 
galaxy (proto-)clusters have the same number of mem- 
bers at all redshifts, a dependence with redshift would be 
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Fig. 7. — Distribution of the densities of IRAC-selected sources 
for the HzRG and RLQ fields in the CARLA sample. A two sided 
K-S test shows no significant difference between the two distribu- 
tions, implying that the environments of the two AGN types are 
similar. 

expected as wc observe galaxies that are about 3 magni- 
tudes fainter than L* at z ^ 1.3 but only galaxies that are 
~ 2 magnitudes fainter than i* at z ~ 3.2. Wylezalek et 
al. (in prep.) uses the CARLA sample to quantitatively 
examine the redshift evolution of the cluster luminosity 
function. The fact that we observe a dependence with 
redshift demonstrates that a non-negligible fraction of 
the IRAC-selected objects in the vicinity of the radio- 
loud AGN are associated with the AGN. If the IRAC- 
selected sources were not associated with the radio-loud 
AGN then their density would likely be independent of 
the AGN redshifts. 

4.6. Dependence on Radio Luminosity 

Figure |9] shows how the density of the CARLA fields 
depends on the radio luminosity of the targeted radio- 
loud AGN. Spearman rank correlation coefficients do not 
show a significant dependence of the CARLA field den- 
sities on LsooMHz- This is in agreement with previous 
studi es at these high radio luminosities (|Donoso et al.l 
|2010[) and is consist ent with expectations f rom the uni- 
fied model for AGN ()Urrv fc Padovani|[l995l) . We discuss 
this result in the context of other studies at lower red- 
shifts in Section 6. 

We also test the difference of the density distributions 
between the HzRG and RLQ fields as a function of radio 
luminosity. We fit a Gaussian to the density distribution 
in four radio luminosity bins (10^^ — 10^* W Hz~^, lO^'^ — 

1028-25 W Rz-\ 1028-25 „ ^028-5 W Hz"!, lO^S'^ - 102^ W 

Hz^^, chosen to provide similarly sized samples) to deter- 
mine the mean and sigma of the histograms. Two-sided 
K-S tests do not show a significant probability for the 
distributions to be different. Spearman rank correlation 
coefficients do not imply a dependence of the density on 
radio luminosity for either of the two types of AGN. 

5. THE ENVIRONMENTS OF HZRGS AND RLQS 

It has long been debated whether the environments 
of radio-quiet AGN are any different from their radio- 
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TABLE 2 
Comparison to previous work on type-1 and type-2 AGN environments. 



Reference 



Sample Size 



(W Hz" 



type-1 vs. type-2 environment 



Yates et al. (_19S9) 
Smitl^^Ieckman 
Falder et al. f2010 
(2010 



Donoso et al 



This work 




25 RGs 


0.15 < z < 0.82 


1 


10 


35 RGs, 31 RLQs 


< 0.3 


5 


10 


75 RLQs, 27 RGs 


z ~ 1 


2 


10 


307 RLQ, ~ 1400 RGs 


0.4< z < 0.8 


2 


10 



166 RGs, ISO RLQs 



1.3 < 2 < 3.2 



3 ■ 10^ 



No difference 

No difference 

Slightly denser type-2 environments 

type-2 more clustered for 10 < I-soomhz < 10 

similar clustering for Lsoomhz > lO" 

No difference 
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Fig. 8. — Surface density of IRAC-selected sources in the HzRG 
(green crosses) and RLQ (orange diamonds) fields as a function of 
redshift. The dark and light grey shaded areas show SgpUDS i Icr 
and EspUDS i Scr density, respectively. A Spearman rank correla- 
tion analysis shows a dependence of surface density with redshift. 



Fig. 9. — Surface density of IRAC-selected sources in the HzRG 
and RLQ fields as a function of radio luminosity. The dark and 
light grey shaded areas show SspUDS i 1"" ^^nd SgpUDS i 3(t den- 
sity, respectively. A Spearman rank correlation analysis shows no 
significant dependence of surface density with radio luminosity. 



loud analogs. Many studies show that radio-quiet 
AGN, which are usually associated with less mas- 
sive spiral galaxies, are less clustered than radio-loud 
AGN (e.g..lHutchings et~aI|[l9 99l:'Ka uffniann et al.ll2008l: 
iHickox et al.ll200a JFalder et al. 201(3). Several theories 
have been proposed to explain this difference between the 
environments of radio-loud and radio-quiet AGN. One 
theory, known as jet confinement, posits that the higher 
IGM densities in regions of higher galaxy densities en- 
hance synchrot ron losses from radio jets , thereby making 
them brighter (jBarthel fc Arnaudlll996D . An alternative 
theory suggests that mergers in regions of higher galaxy 
density will produce more rapidly spinning supermassive 
black holes, which will then be more capable of power- 
ing a radio jet (e.g.. lWilson fc Colbertl[l995l : iSikora et alj 
|2007[ ). Regardless, observational results suggest that 
radio-quiet AGN reside in less dense environments while 
the difference between RLQ and RG environments still 
remains an open issue and observations show inconclu- 
sive results. 

From a study of ^ 30 r adio-selected AGN at z < 0.3, 
iSmith fc HeckmanI ()1990[ ) finds that RLQs and RGs at 
z < 0.3 were simila r ly clu stered, consistent with the re- 
sults of lYates et al.l (|1989D based on a comparably sized 



sampl e at similar redshift (see Tab. [2]). iFalder et al.l 
()2010D studies the environments of 75 RLQs and 27 RGs 
at z '^ 1, and finds a small, but notable trend of den- 
sity with AGN type, with the RG environments being 
slightly denser than the RLQ environments. He sug- 
gests, howeve r, that this ni i ght be due to the different 
sample sizes. iDonoso et al.l ()2010[ ) computes the cross- 
correlation signal between a sample of ^ 14,000 radio- 
detected AGN and a sample of 1.2 million luminous red 
galaxies in the SDSS. They also study a sample of r^ 300 
RLQs cross-correlated with the same sample of luminous 
red galaxies. When c omparing the two clustering signals, 
IDonoso et al.l ()2010D finds the RGs are more clustered 
than the RLQs for lO^^ < Li.4ghz < 10^*^ W Hz"!, but 
clustered in a similar way for Li.4ghz > 10^^ W Hz~^. 

We caution the reader that in each study different 
lower limits and definitions of radio-loudness are being 
used. In Tab. [2] we summarize the results of these works 
and list the redshift ranges and the lower limit on radio 
luminosity used to define a radio- loud object. The radio 
luminosities have been converted t o J^boomhz assuming 
a rad io spectral index of a = —0.7 (jMilev fc De Breuca 
I2008D . At lower radio luminosity, Lsoomhz < 10^^ W 
Hz~-^, the studies show somewhat conflicting results. A 
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concerning problem is the small sample sizes providing 
poor statistics and with imperfectly matched samples. In 
iFalder et alj (I201CID . for example, the RGs were slightly 
more radio luminous than the RLQs. 

At low radio luminosities, our understanding of the en- 
vironments of RLQs and RGs is not yet clear. At higher 
radio luminosities, isooMHz > 10^^ W Hz~^, previous 
studies and our analysis agree that there are no signifi- 
cant differences between the environments of type-1 and 
type-2 AGN. These results are consistent with the popu- 
lar orientation-driven (e.g., obscuring torus) AGN unifi- 
cation scenario for the most powerful type-1 and type-2 
radio-loud AGN. However, at lower radio luminosities, 
other factors, such as AGN feeding mechanisms, appear 
to come into play causing less powerful radio- loud type-1 
and typc-2 AGN to, on average, reside in different envi- 
ronments. 

6. SUMMARY 

We have studied the environments of 387 radio-loud 
quasars and radio galaxies at 1.3 < z < 3.2 that have 
been observed at 3.6/im and 4.5/im with the Spitzer Space 
Telescope. Our study clearly confirms that radio-loud 
AGN reside in rich environments. By applying a color 
cut, we identified likely high-redshift galaxies, which we 
refer to as 'IRAC-selected sources', and computed their 
surface density in a 1' radius around the radio-loud AGN. 
More than half of the fields, 55.3%, are overdense at the 
2 a level when compared to the blank field density of 
IRAC-selected sources derived from the SpUDS survey. 
If we instead consider a 5 ct overdensity, 9.6% of the 
fields are overdense. We detect a significant rise in sur- 
face density of IRAC-selected sources towards the po- 
sition of the radio-loud AGN in contrast to a homoge- 



neous radial surface density distribution for the SpUDS 
apertures. This reassures us of an association of the ma- 
jority of the IRAC-selected sources with the radio-loud 
AGN. These CARLA results provide a promising list of 
high-redshift galaxy (proto-)clustcr candidates for fur- 
ther study. 

We also study the dependence of environment on ra- 
dio luminosity. No statistically significant correlation is 
observed. We do, however, find a significant dependence 
on redshift, with higher redshift fields being slightly less 
dense. This result is consistent with galaxy cluster evo- 
lution and reassures us that most of the IRAC-selected 
sources are indeed associated with the targeted radio- 
loud AGN as otherwise we would have been unlikely to 
have observed a redshift dependence. In Wylezalek et al. 
(in prep.) we discuss how this trend can be explained 
in the context of the evolution of the luminosity/mass 
function of galaxy clusters. 

We do not find a significant difference between the 
RLQ and HzRG fields, neither in their overall distribu- 
tion nor in their dependence with redshift or radio lu- 
minosity. We conclude that for high radio luminosities, 
-^sooMHz > 10^^-^ WHz~^, orientation-driven unification 
models appear valid. 
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